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Abstract 

A quantum dot driven by two ac gate potentials oscillating with a phase lag may be regarded 
as a quantum engine, where energy is transported and dissipated in the form of heat. In this 
chapter we introduce a microscopic model for a quantum pump and analyze the fundamental 
principle for the conservation of the charge and energy in this device. We also present the basics 
of two well established many-body techniques to treat quantum transport in harmonically time- 
dependent systems. We discuss the different operating modes of this quantum engine, including 
the mechanism of heat generation. Finally, we establish the principles of quantum refrigeration 
within the weak driving regime. We also show that it is possible to achieve a regime where part 
of the work done by some of the ac fields can be coherently transported and can be used by the 
other driving voltages. 
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I. INTRODUCTION 



In 1824 the french physicist Nicolas Leonard Sadi Carnot, better known as Sadi Carnot, 
in "Reflections on the Motive Power of Fire" settled the principles of the modern theory 
of Thermodynamics by pointing out that motive power (concept later identified as work) is 
due to the fall of caloric (concept later identified as heat) from a hot to cold body (working 
substance). These ideas, that provided the scientific support for the technological jump 
based in the steam engine were not well understood at that time. They were actually 
discovered and further elaborated thirty years later by the German Rudolf Clausius and the 
British William Thomson (Lord Kelvin). The fundamental principles ruling the operation 
of the thermal engines were then summarized in the two basic laws of Thermodynamics. 
While the first law simply stresses the conservation of the energy, the second one deals with 
the subtle distinction between a kind of energy that can be used and another one that is 
dissipated in a physical process, as well as on the balance between both of them. 

In 1851 Lord Kelvin also discovered Thomson effect, and showed that it was related to 
other thermoelectric phenomena: Peltier and Seebeck effects. Unlike Carnot 's machines, 
these effects are related to non-equilibrium processes. However, they also bring about the 
conceptual distinction between different kinds of energies: one that is transported in some 
direction due to a voltage or temperature gradient, but being different from the Joule heating 
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in the sense that the first one is reversible while the latter involves dissipation and, thus, 
irreversible effects. 

Nowadays, we are witnessing a technological trend towards an increasing miniaturization 
of the electronic components. This is accompanied by a significant activity within com- 
munities of the basic sciences, in the search of a better understanding of the behavior of 
materials and devices with sizes in the range between Inm to 10/im as fundamental pieces 
of electronic circuits. Paradigmatic examples are the quantum dots fabricated in the inter- 
faces of semiconductor structures where confining gates for electrons and circuits are printed 
by means of nanolitography within an area of a few /im^ (see Fig. [T]). Due to the small 
scale of these systems, they present some physical features that resemble the molecules. In 
particular, the landscape of their spectra contains well defined quantum levels where elec- 
trons propagate almost perfectly preserving the phase of their wave functions. However, 
they are not isolated from the external world but coupled to the substrate, gates, wires and 
external fields that induce the transport of electrons. For this reason, they are classified as 
"open quantum systems" that operate out-of-equilibrium conditions. The "external world" , 
instead, contains pieces that act as macroscopic reservoirs with which the "small quantum 
systems" exchange particles and energy. Due to the mixed nature of these systems, the con- 
cepts of classical electrodynamics and thermodynamics cannot be simply applied to them 
and theoretical tools that are amenable to capture their quantum properties as well as the 
coupling to their environment are necessary. 

In this chapter we focus in a particular kind of devices named "quantum pumps" . They 
have been realized experimentally, precisely, in quantum dots, where ac voltages are applied 
at their walls, and a current with a dc component is generated in the absence of a stationary 
voltage difference (see Figure 1). From the theoretical point of view one of the interests 
in these systems is that they can be described in terms of simple Hamiltonians, including 
macroscopic pieces that represent the wires connected to the quantum dot as well as explicit 
terms for the time dependent forces that induce the transport process. The latter feature 
makes an important difference, in comparison with setups where the transport is induced 
by means of a stationary voltage difference. This is because in such systems the work of 
the forces that keep such a bias fixed is not explicitly taken into account in the model 
Hamiltonian but introduced as a boundary condition. 

We can imagine situations in which the macroscopic wires connected at the quantum dot 
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FIG. 1: Quantum dots fabricated with nanolitography in a semiconductor interface. The two leads 
going upwards and downwards make the contact between the dot and external reservoirs. The 
gate at the left can be used to apply a voltage to introduce a rigid shift of the positions of the 
energy levels of the dot. The two extra leads connected at the right hand wall of the dot can be 
used to apply ac potentials at those points. Figure by Charles Marcus, Mesoscopic Lab, Harvard 
University. 

are at different temperatures. Additionally, the time dependent forces do make work on 
the system. Thus, we can regard the quantum pump as a microscopic engine where heat 
can be exchanged between two sources at different temperatures, while work is provided 
to the system and part of the energy is dissipated. As we will show, such an engine could 
even operate as a refrigerator, where there is a net heat flow from the reservoir at the lowest 
temperature to the one at the highest one. The goal of this chapter is to introduce theoretical 
tools for the analysis of the fundamental conservation laws at the microscopic scale, the 
explicit evaluation of the power developed by the intervening forces and the distinction 
between reversible flows and dissipated energy. 

This chapter is organized as follows. In the next section we introduce a simple microscopic 
model for a quantum pump and we will discuss the fundamental principle of the conservation 
of the charge and the energy in this device. In section III we introduce the basics of two 
well established many-body techniques to treat quantum transport in harmonically time- 
dependent systems. The first one starts from the explicit microscopic Hamiltonian for the 
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system, forces and environment which is solved by recourse to non-equihbrium Green's 
functions. While the second one is based in the notion of scattering processes that the 
electrons experience as they cross through a quantum system under ac driving. Although 
this chapter is self-contained, we will not include a complete tutorial on these two techniques 
but we adopt a practical point of view, presenting just the main ideas while we defer the 
reader to more specific literature on many-body techniques for further details. In section 
IV we present explicit expressions for the energy flows in terms of the Green's functions 
and the scattering matrix elements introduced in section III. We also discus the nature of 
the different components contributing to the total energy flow. Section V is devoted to a 
summary of the different operating modes that we were able to identify in our quantum 
engine. Finally, in section VI we conclude with a discussion of the possible directions to 
extend these ideas. 

II. BACKGROUND 
A. Model 

We start by defining explicit Hamiltonians to describe the quantum electronic system 
as well as its environment. For the sake of simplicity, let us focus in a system like the one 
sketched in Fig [2] where the central system C, on which the time-dependent forces are acting, 
is placed between two wires: one located at the left (L) and the other at the right (i?) of C. 
We can identify C with the quantum dot of Fig. [1], with the two ac voltages applied at the 
two extra leads connected at the wall. The L and R reservoirs of Fig. [3] correspond to the 
up and down ones in Fig. [H 

Along this chapter, we make the following simplifying assumptions on the system: (i) We 
shall take into account only the two wires as the external environment for the central system 
and disregard other effects like, for example, the influence of the phonons of the substrate. 
Such a simplification is expected to be reasonable if we concentrate on the behavior at 
sufficiently low temperatures, (ii) We also assume that the electrons do not experience 
any kind of many-body interactions, like the Coulomb repulsion between electrons. This 
assumption is justified in the description of the metallic wires where we can expect an 
efficient screening but it is justified within the system C only when the structure is strongly 
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FIG. 2: Sketch of the setup. The central system is connected two two infinite wires, which play 
the role of macroscopic reservoirs. In this example, the central system contains a profile of two 
barriers at which ac fields that oscillate with the same amplitude Vq and frequency l^o, and a phase 
lag 62- di. 

connected to the wires and allows for the screening of the wires to penetrate into it. In this 
context, the spin degrees of freedom of the electrons behave independently of one another. 
For this reason, in order to simplify the notation, we do not consider them explicitly. In the 
case of considering them we must simply write a factor 2 in front of the final expressions for 
the currents and densities. The full Hamiltonian reads: 

H{t) = Hl + Hc{t) + Hn + H,, (1) 

where 

Ha = ^ ^ ^fcaC^gCfcQ, , (2) 
ka 

being a = L,R, are Hamiltonians of free spinless electrons that represent the wires. We 
stress that these systems are macroscopic, i.e. they contain a very large number of degrees 
of freedom and are in thermodynamic equilibrium. This means that they are completely 
characterized by their density of states 

Paiuj) =271^^5 (U - £ka/fi) , (3) 
ka 
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and the Fermi distribution function: 

/-(^) = efe(^-^«) + 1 ' 

with 13a = ksTa. The second term, Hc{t). describes the quantum structure under consider- 
ation as well as the time-dependent gate potentials acting on it. The ensuing Hamiltonian 
depends on the geometry of the structure as well as on the interactions that we want to take 
into account. In the absence of many-body interactions, this system may be described by a 
single-particle Hamiltonian of the form: 

Heir, t) = + U{r) + <5(r - R,) eVi{t) , (5) 

^ 1=1 

being m the mass of an electron, which corresponds to a finite number of time-dependent 
potentials that we assume have the simple single-harmonic dependence: Vi{t) = Vq cos(r2o^ + 
6i). The potential U (r) contains the information of the confining walls, barriers and defects of 
the structure. For a detailed discussion of the conservation laws and for treating the problem 
with the Green's function formalism, it is convenient to express this Hamiltonian in second 
quantization. To this end we must define an appropriate single particle basis to represent 
the relevant operators. As the structure under study occupies a reduced region of the space, 
it is comfortable to work with a single particle basis that is labeled by spacial coordinates, 
like that defined by the Wannier functions. It is, thus, useful to work on a discrete lattice 
containing a finite number (A^) of sites and a basis of single-electron states that are localized 
on the lattice positions. The resulting Hamiltonian in second quantization corresponds to a 
tight- binding model. For simplicity, we consider this model in one-dimension (ID), although 
this is not an essential assumption: 

N N-1 

Hc{t) = YJ\.^i + eVi{t)]c\ci - w 5^(cJq+i + H.c) , (6) 
1=1 1=1 

where the term with w = {l\{fi^ /2m)d'^ / dx'^\l + 1), being |/) single-electron basis state 
localized at the lattice position describes the kinetic energy of the electrons through 
jumping processes between nearest-neighbor sites of the underlying lattice. The term with 
Si — {l\U{x)\l), defines a static energy profile for the structure: it contains the information of 
the existence of barriers and wells. For a system with impurities, this profile can be defined 
in terms of a random amplitude and this model reduces to the Anderson model. The term 

8 



with Vi(t) = {l\ J2j=i^{^i ~ ^j)^o cos{Qot + Sj)\l) represents the time dependent gates, being 
finite at the M pumping centers and vanishing otherwise. Finally, the term He describes 
the contacts between the central system and the reservoirs. In our simple ID model for the 
central structure the L lead is connected to the first site, 1 = 1, and the R to the last site, 
I = N, of the central structure. The Hamiltonian reads: 

He = -WeL ^{cl^Ci + H.C) - WeR ^(4rCjv + H .c) , (7) 
kL kR 

which describes hopping processes between the states ka within the wires and the points of 
C at which the contact between the two systems is established. 

Before closing this subsection, let us mention that, depending on the physical problem 
under consideration, there may be other time-dependent terms in the Hamiltonian. For 
systems with ac voltages applied at the L and R wires, we should consider a dispersion 
relation with a time- dependent component in addition to the static one, . By recourse 
to a gauge transformation, it can be seen that this type of ac voltages can, equivalently, 
be included in a time-dependent phase in the contact hopping Wc (see Jauho et al. 1994). 
Another possible time-dependent term is that originated by an electric field derived from a 
time-dependent vector potential A(t) (see Arrachea 2002). That physical situation would 
take place, for example, when the central system is bended and closed into an annular 
geometry threaded by a time-dependent magnetic fiux. In terms of the Hamiltonian this 
introduces a shift: p p — (e/c)A(t) in the momentum of the electrons, being c the velocity 
of light. In the tight-binding basis this translates into time-dependent phases in the hopping 
parameter w of Hc{t) along with the periodic boundary condition + 1 = 1. Finally, 
the study of the coupling to external classical radiation fields is usually treated within the 
so called dipolar approximation, which in our second quantization language results in a 
diagonal voltage profile Vj{t) as in Eq. ([6]) with Vj ~ jVo, with Vq constant and 6j = 6,Wj 
(see Kohler et al. 2005). 

B. Conservation laws and instantaneous currents 

1. Particle currents and the conservation of the charge 

A consistent way to define expressions for the electronic currents along the different 
pieces of the structure is starting from the evolution of the electronic density rii = cjci. The 
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variation of rii is due to the difference between the charge flow exiting and entering the 
infinitesimal volume that encloses that point: 

- ej^{ni) = Ji{t) - Ji^i{t) , (8) 

where Ji{t) denotes the current exiting the site I towards the neighboring site / + 1. We 
denote with a positive sign the flows pointing from left to right. The variation in time of the 
local charge can be calculated within the Heisenberg picture by recourse to the Eherenfest 
theorem: 

ej^{n,) = --^{[H,c\(t)cm)- (9) 

Thus, the explicit evaluation of the above commutator defines an explicit expression for the 
current. If we consider a site within C we obtain: 

m = ^(cKt)Q+i(t) - ciAtym . (lo) 

It is easy to verify that the above expression coincides with the mean value of the operator 
ev, being v the velocity expressed in second quantization in the basis of localized functions. 
Similarly, if we consider the contact between C and one of the reservoirs we get the following 
expression for the current that exits the reservoir a: 

Ut) = ^ W - clit)cUt)) , (11) 



h 

ka 



where la = 1,N for a = L,R. 



2. Energy currents, power, and the conservation of the energy 

In order to define energy currents, we proceed along a similar line as before. In this 
case, we analyze the evolution of the energy density at a given elementary volume and write 
the equation of the conservation of the energy. As our Hamiltonian Hc{t) contains terms 
involving positions up to nearest-neighbors, the smallest volume for analyzing the evolution 
of the energy density in our ID lattice is that enclosed by a box confining a bond of nearest- 
neighbor sites (see Figure [3]) • If we denote by the total energy stored within such a 
box, the equation for the conservation of the energy reads: 

= JFAi) - Jfit) + m) + ^m(t) , (12) 
10 




FIG. 3: Analysis of the energy balance in our setup. The dashed boxes enclose the elementary 
volume of the system where the evolution of the energy is studied. It contains two sites of the 
underlying lattice. The arrows indicate the direction that we have defined as positive for energy 
currents along the different pieces of the system as well as the powers done by the fields. 

where the first two terms denotes the difference between outgoing (from / + 1 to / + 2) and 
incoming (from / — 1 to /) flows, with the same sign convention as in the case of charge flows, 
while the last two terms denote the power done by the external fields, which are defined as 
positive when it is provided by the forces. The latter terms vanish if the time-dependent 
gate potentials are not acting at the points / and / + 1 enclosed by the box. Our box can 
also enclose the contact bond between the reservoir a and the central system (see Fig. [3]), 
in which case we get: 



dt 

Ji{t)-J^_,{t) + P^{t), (13) 



dt 

where J^it) denotes the energy flow that exits the reservoir a. As in the previous subsection, 
the explicit expressions for the flows and powers are derived by recourse to the Eherenfest 
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theorem: 

= -^([■ff.(':i(*)cK*)<;iW + ';i+iW<;;+iW<:i+i(*)-«'{cK*)Q+iW + <;|+i(f)qW})]> 

+e^{c|(t)c,(t)> + (4,(t)cH,(t)> , (14) 

with ei{t) — s'l + e.Vi{t), and a similar expression for the volume enclosing the contact bonds 
between C and the reservoirs. Prom the evaluation of the previous terms, we obtain the 
explicit expressions for the energy currents: 

Jiliit) = y [^^(cJ(t)Q+2(t) - cl,{t)ci{t)) - ei+,{t){cl,{t)ci+2{t) - cl,{t)cl+^m , (15) 

Ja(t) = '-^J2'ka{cl{t)cia{t) - cL(t)c,.(i)) , (16) 
ka 

and power developed by the ac voltages: 

m) = e^iclitMt)) . (17) 



C. Continuity equations, dc charge and energy currents and mean powers 

In the absence of sinks and sources, the average of the charge enclosed by any volume of 
the sample over one cycle with of period r = 27r/ilo must remain constant, which defines 
the following continuity equation for the dc charge current (microscopic Kirchoff law) : 

Ji^Ji'^J, (18) 

for arbitrary /,/' along C, where we denote A = 1/r dtA{t). 

Analogously, for any volume enclosing lattice points running from / + 1, we get: 

I' 

jf = J?+J2P^^ (19) 

where the last term defines the power done by all the voltages enclosed by the volume. The 
above equation reduces to a "Kirchoff law" for the dc energy current when we enclose a 
region that is free from the time-dependent voltages, in which case the last term vanishes. 

D. Heat current 

— E 

At this point, it is important to mention that the dc energy current J defined above 
does not necessarily coincide with the heat current which we will denote J^. This is because 
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what is understood as "heat" is usually the energy transferred from one system to another 
as a consequence of a temperature difference. In order to understand this difference, let us 
consider that our reservoirs are at temperature T = 0, let us place our volume enclosing 
the contact between one reservoir and the central system and let us assume that we are not 
enclosing an ac local voltage within it. Now, let us analyze the following picture based on 
heuristic arguments that we shall better formalize it in subsequent sections. Let us assume 
that the ac voltages are so weak in amplitude and oscillate with such a low frequency that 
we can disregard the power they develop. Let us assume that, anyway, they are able to move 
a small portion of electrons with energies very close to the Fermi energy of the reservoirs /i, 
that we recall is the same for L and R. This weak motion give place to currents of particles 
and energy which, from the definitions flTT]) and flTBl) . are approximately related through 
J^{t) ~ {fi/e)Ja{t). The same relation holds for L as well as R reservoirs and a small 
current of particles may translate in a large current of energy, since fi can be large. The 
dc component of the charge current may be finite and should be positive in one reservoir 
and negative in the other one, indicating that there is a net flow of charge between L and 
R or vice versa. The above relation, therefore, tells us that there is a concomitant net flow 
of energy from a reservoir to the other one, in spite of the fact that we are assuming that 
both reservoirs are at zero temperature. One could complain that we have not taken into 
account the power done by the time-dependent fields which would tend to heat the system. 
However, let us recall that an appropriate choice of fi would allow us an arbitrary large 
value of the energy flow, against which we can disregard a contribution like ( fTTl) for weak 
and slow ac voltages. In summary, the above considerations lead us to conclude that there 
may exist a net energy flow which cannot be identified as "heat" but has rather a convective 
nature. In order to better quantify heat, it is thus natural to subtract from the energy flow, 
the convective component (/i/e)J, i.e.: 

Jf^ = jf-^^J, 4« = jf-^J. (20) 

Notice that, while the convective term is constant along all the pieces of the system due to 
the continuity of the charge, the heat and energy currents may have different values due to 
the contribution of the power done by the external voltages, see Eq. (1191) . 

To give a more formal definition of the heat flow , let us consider the particle and 
the energy balances for a given reservoir a that is kept at fixed both the chemical potential 
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Ha and the temperature T^. If the charge current J a and the energy flow enter the 
reservoir a, then the number of particles and the energy of a reservoir should change. This, 
in turn, would change the chemical potential and the temperature of a reservoir. In order 
to keep Ha fixed, the electrons have to be removed with the rate Ja/e out of the reservoir. 
Therefore, while keeping constant we necessarily remove energy with rate ^aJa/^- We 
stress that the convective energy HaJa/^ is taken out at equilibrium conditions, therefore, it 

— E 

can be reversibly given back. In general this energy does not coincide with energy flow 
entering the reservoir. To prevent heating of reservoir one needs additionally to take out the 
energy with rate = ~ ^J^aJ al(^ without taking out particles. Since the reservoir can not 
produce work, the only way to remove the remaining energy is to put it in contact with other 
large body playing the role of a thermostat. The energy exchange between the reservoir a 
and the thermostat is essentially irreversible. For this reason we interpret this part of the 
total energy as "heat" and we identify as the heat flow. If the thermostat would be absent 
the temperature of a reservoir would change. As we will show can be directed either to 
the reservoir from the central system or back. Hence, in the absence of a thermostat, the 
reservoir can be either heated or cooled. We stress that the energy transported at the rate 
becomes "heat" only deep inside the macroscopic reservoir, where the electrons scattered 
by the dynamical central system, are able to equilibrate. 

III. STATE OF THE ART 

In this section we briefly review the many-body techniques to evaluate the currents and 
powers defined in the previous section. 

A. Green's functions formalism 

1. Expectation values of observables and Green's functions 

The expectation value of any one-body observable, {A{t)) — '^ii'{l'\A{t)\l){cl,{t)ci{t)), 
can be regarded as follows: 




(21) 



1,1' 
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being 



G< 



.(t,t')^i{cl(t')ci(t)), 



(22) 



a "lesser'' Green's function. Our goal, now, is to derive equations for the evolution of this 
Green's function and strategies to solve them. 

2. Brief review of the theory of the non- equilibrium Green's functions 

The formal theory of non-equilibrium Green's function has been developed independently 
by Kadanoff and Baym [Kadanoff and Baym 1959], Schwinger [Schwinger 1961] and Keldysh 
[Keldysh 1962]. The structure of that theory is very similar to the one of causal Green's 
functions at zero temperature (see Mahan 1990), except for the fact that in non-equilibrium 
situations, the assumption that the state of the system at time -|-oo differs just in a phase 
from the state in — oo does not longer holds. The way to overcome this inconvenience 
is to define the evolution along a special contour C that defines a round trip, first going 
from — oo to -|-oo and then going back to — oo. As in equilibrium problems, the precise 
description of that evolution can be accomplished with the help of Wick's theorem and 
Feynman diagrams and one of the big powers of this technique is the possibility of treating 
many-body interactions in a systematic way. 

We skip here the technical just highhghting the main ideas leading to some useful identi- 
ties, and we defer the reader to more specialized hterature (see Mahan 1990, Haug and Jauho 
1996, Caroh et al. 1971, Pastawski 1992, Jauho et al. 1994). Instead of the time-ordering op- 
erator used in the equilibrium theory, it is convenient to work with contour-ordered Green's 
functions: 



where 1,1' is a schematic notation that labels the electronic degrees of freedom and time 
in the same index and the operator Tc denotes time-ordering along a contour that begins 
in — oo evolves to -|-oo (C+) and then turns back from +oo to — oo (C_). This function 
corresponds to the casual "time-ordered" , the "lesser" , the "anti-time-ordered" or "greater" 
function depending on the position of the two times along the closed time-contour: 



^(1, 1') = G,{1, 1') , ii, ty e C+ , ^(1, 1') = G<{1, 1') , ii, e C+ , ty e C_ , 
G{1, 1') = Gc{l, 1') , ti,tr e C_ , Gil, 1') = G>il, 1') , t,, e C^M' e C+ , (24) 



G(l,l') = -^(^c[c(l)ct(l')]), 



(23) 
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with: 



1') = -I &{ti - tr){ci{h)clit'i)) + ^ Q{ti' - h){c\{t[)c,{U)) , 
G^{1, 1') = -z e(ti. - ti)(ci(ti)4(t;)) + t Q{h - tr){c\{t[)ci{h)) , 

G<(i, 1') = ^ (4(t;)ci(to) , G>ii, 1') = -I (ci(to4(t;)) , (25) 

which are not independent functions, but satisfy: Gc + Gc = + . It is also convenient 
to define "retarded" and "advanced" functions: 

= e{h-tr)[G>il,V)-G<{l,l')], 
G^(l,l') = e(ti,-ti)[G<(l,lO-G>(l,l')], (26) 

which are also related through G^ -G^ = G> -G<, G^{1,V) = 1)]* and G'<(1, 1') = 

-[G>{i',i)r. 

In our simple model of non-interacting electrons, we can split by convenience the Hamil- 
tonian in two parts H{t) = HQ{t) + H'{t), with both terms being of one-body type but Hq 
being easily solved. The evolution of this Green's function is given by Dyson's equation, 
which for a one-body Hamiltonian reads: 

^(1, 1') = Go(l, 1') + h-^ d^X2dt2G{l, 2)H'{2)Go{2, l') . (27) 

Notice that the above equation actually represents a matricial integral equation if we distin- 
guish the positions of the times as in fl25l) . A convenient tool to derive explicit equations for 
the different components (125|) is a theorem due to Langreth, which states given a product 
of contour-ordered Green's functions of the form: 

Gih,h,) = ^ dt2Giih, h)G2it2, ty) , (28) 

then, the following relations hold for the different components: 

G^'^(tl, tv) = r dt2Gf'^{h, t2)G^'^{t2, tv) , 

/— oo 
dt2[Gf{ti, t2)G<'>{t2, ty) + Gf •>(ti, t2)G^(t2, ty)] . (29) 
-oo 

Therefore, if we want to compute G^, in order to evaluate expectation values of observables, 
we must also solve two coupled equations for that function and G^. 
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3. Green's functions and Dyson's equations in our problem 



Let us first split our Hamiltonian as follows: -ffo(^) = Hl + Hn + Hc{t) and H' = He- 
The Dyson's equation for the retarded function reads: 

«,,.(MO^.L^(M',.r.i:/*«„.(Mo..L.(.,o, (30, 

where gj ,j,{t,t') is the contour-ordered Green's function of Hoit) and run over all the 
electronic degrees of freedom of this Hamiltonian. If we write the above equation explicitly 
for one of indexes in the reservoir and / G C: 

Gi,kait, t') = -Wc ^ dt^Gi^Ut, t,)gl^^,^{t,,t') , (31) 

For the two indexes 1,1' & C: 

Gi,v{t,t') = gl,{t,t')-w, J2 f dtiGi,ka{t,h)gl^,{t,,t') , 

/oo 
dtidt2Gi^UtM)^a{tiM)9lAt2,^') ' (32) 
_ -oo 



a=L,R- 



where going from the first to the second identity we have substituted fl3T|) . We have also 
defined the "self-energy" : 

S«(tl,t2) = kca|'$^<,fcjtl,t2) . (33) 
ka 

In order to evaluate the currents that we have defined in section (III Cp . we need Gfii{t, t') 
for 1,1' E G and G'^^i^(t,t'). Applying the Langreth's rules fl2^ to the above equations we 
get (see Arrachea 2002 and 2005): 

/ + 00 
dhlG^iait, ti)9'uX^h,t') + G<f (t, ti)^7i:i.Jti, t')] , (34) 
-oo 

along the contact and 

/ + 00 
dt,dhGf^,^(t, h)^<'>{h,t,)GiAt,, t') , 

(35) 

for coordinates /, /' belonging to the central system G. The latter equation is obtained 
after some algebra and after dropping a term that contains g^'jf{t,t') which can shown to 
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relevant only in the description of transient behavior (see Jauho et al. 1994). The different 
components of G^{t,t') within the reservoirs are straightforwardly evaluated: 

alikaituh) = -iQ{ti -t2)exp{-i{eka/fi){ti -t2)} , 

9kaMitut2) = A<'>(£fc„)exp{-^(5,jn)(ti-t2)}, (36) 

with Xai^ka) = ifai^ka) sjid X^i^ka) = —i [I — fai^ka)]- With these functions, it is possible 
to obtain expressions for the different components of (133|1 in terms of the density of states 
of the reservoir Pa(co') given in ([3]): 

/+00 7 
— e-^-(*-*^)r,(o.) , 
-oo 

/+00 7 
^ |^e-(*-*^)A<->(;k^)r„(u;), (37) 

with Ta{uj) = \Wca\'^Pa{uj). 

To evaluate flMl) and fl35l) we still have to calculate the retarded function within the 
system C (recall that the advanced function can be obtained from Gfi,(t,t') = [G'/?^(t', t)]*). 
The equation for the retarded function is the retarded component of fl52]) and can be derived 
by applying Langreth rules on this equation. The result is: 

G5,(t,t') = glf{t,t') + h-' rrfiirfi2GjJt,ti)Sf (ti,t2)C'(*2,0 , (38) 

with g^'i^{t,t') being the retarded Green's function of the system described only by Hc(t) 
isolated from the reservoirs. This function is, in turn still an unknown in our problem. 
Instead of explicitly evaluating it, we find it more convenient to operate with (1381) in order 
to find an equivalent equation for G^{t, t') as follows. We first derive an equation of motion 
for g^'i^it, t'), starting from the very definition of the retarded function, see Eq. (1261) . and by 
writing the evolution of Cj{t) with Hc{t) in the Heisenberg representation: 

zhcj (t) = [Hc{t),cj{t)], (39) 

we get: 

- th^g^'^t, t') - g^'^t, t')Hc{t') = hU{t - t') , (40) 

where g^'^(t,t') is a. N x N matrix with elements g^\^(t,t') and i is the N x N identity 
matrix. The above equation means that {-ihd/df -Hc{t')} = [g^'^]'\ Therefore, we have 
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not evaluated explicitly the function g^'^{t,t') but we have identified an operator which is 
its inverse. We act with this operator from the right of (|38|) and we consider the following 
splitting of the central Hamiltonian Hcit) = + H'^(t), where Hcit) collects all the 
explicit time-dependent terms of Hc(t) and Hq, the remaining ones. We get: 

-ih-^G'^{t,t')-G'^{t,t')H^- J dtiG'^{t,ti)f:'^{ti,t')-G'^{t,t')H'c{t') = hl6{t-t') , (41) 

being T,ji^i,(t,t') = J^a^iM^i'M^ai't^t')^ the matrix elements of T,{t,t'). We define a function 
G°'^{t,t'), such that = {-ihd/df - - ±^}: 

- ih-^G^'''{t, f) - G'^'^^it, f)H^c - j dtiG''''^{t, ti)S^(ti, t') = hi6{t - t') . (42) 

Multiplying (155]) by the right with G^'^, we finally find the following equation for the full 
retarded Green's function 

^^(t, t') = G'^^^t, t') + h-^ f dt^G'^it, t^)H'c{ti)G°''^{h, t') . (43) 

Jt' 

This equation is completely equivalent to 0381) but has the advantage that the function G^'^ 
is an equilibrium Green's function, which evolves according to the stationary terms of the 
full Hamiltonian H . The above expression has a structure which is particularly adequate for 
perturbative solutions in the time dependent part of Hc{t). We shall exploit this property 
later. 

Equation (l42l) can be easily solved by performing the Fourier transform: G^'^{u) = 
h-^ £ dTe-'^^+'''>G'^'^{T), with r] > 0, since, as we have mentioned before, it corresponds 
to an equilibrium Green's function that depends on t — t'. The result is : 

G°'^(cu) = [h{uj + i7])i -H^c- S'^(^)]~' , (44) 

and can be explicitly evaluated by simply inverting the above N x N complex matrix. 
Substituting in fH5]) and performing a Fourier transform in t — t', Eq. fH5]) results for our 
specific Hamiltonian ([6]): 

G'^it, oo) = G°'^(cu) + J2 e-'''''"'G'^{t, oo + kno)V{k)G^''^{uj) , (45) 

k=± 

where the matrix V{1) contains elements Vi^i' = Si^i' J2jLi^i^i " Xj)eVoe~^^^ and V{—1) = 
[V{1)]*. The linear set (H5l) has the same structure as the dynamics of a problem in which 
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electrons with a given energy hoo interact with a potential V emitting or absorbing an energy 
quantum hQo and scatter with a final energy huj ± hQo. The solution of fH5|) leads to the 
complete solution of the problem. 

Due to the harmonic dependence on the time t of these equations, the retarded Green's 
function can be expanded in a Fourier series as follows: 

+00 

G^(t,cu) = ^ ^(n,cj)e-^"^«*. (46) 

n=— 00 

We give the name of Floquet component to the functions Q{n, cu), because fl46p has a similar 
structure as that proposed by Floquet for the wave functions of time-periodic Hamiltonians. 
The different components obey the following useful identity: 

Q{n, U) - G\-n, LUn) = ~* X] ^(^ + '^'j UJ _n')t {uJ -n')Q {n' , CU-n')^ , (47) 

n' 

where we have introduced the following notation Un = uj + nQo. To prove fH7|) we start from 
the definition (!26l) of the retarded Green's function for indexes /, /' G C. Replacing (!35|) and 
inserting there the representation ( H6i) . we get: 



^^-iHt-t')+no{k,t-k,t')]g^f^^^ ^^Y{tu)g\k2, u) , (48) 

where T{uj) contains as matrix element ri^i/{uj) = Si^iaSi'^iaTai^)- Calculating the Fourier 
transform of this function with respect to t — t' and collecting the n-th Fourier coefficient 
()46|) we find: 

g(„..)^Vrj^g(" + "'--')rWK^')'. ,>o. (49) 



which leads to the identity fl47j) using: 



^ V I } - tTT6{uj ~ 00') . (50) 



00 — 00' + trj [uj — 00' 

It is important to remark that, for Vq = the identity fj47|) reduces to the following identity 
between equilibrium Green's functions: 

p{lu) = G'°'^(^) - G°'^^(^) = -zG'°'^(^)f(^)G'°'^(^)t . (51) 
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4- Perturbative solution of the Dyson's equation 

Sometimes, in order to derive analytical expressions, it is convenient to solve the set (HSj) 
by recourse to a perturbative expansion in V (see Arrachea 2005). The solution up to second 
order in this parameter is obtained by writing fl45p evaluated at Un, for n = —2, . . . , 2 and 
back-substituting the equation evaluated at uj2 into the one evaluated at ui, and the latter 
into the one evaluated at uj, and a similar procedure with uj_2 ~^ ^-i and the latter into uj. 
If we then collect all the coefficients of e~*"^o* in the resulting expression and recalling the 
representation (j^6|) . we obtain 

+2 

G^{t,uj) ~ J2 ^(^,^)e""'^°*, (52) 

n=-2 

with: 

k=±l 

g{±2,uj) = G°'^(cu±2)V(±l)^(±l,u;). (53) 
The reader can easily extend the procedure to evaluate higher order terms. 

B. Scattering matrix formalism 

To calculate the charge and energy flows generated by the driven central system in the 
wires one can also use the scattering approach. Within this approach we consider the central 
system as some scatterer which reflects or transmits electrons incoming from the wires. The 
electrons coming, for instance, from the left wire can be transmitted to the right wire or 
can be reflected back to the left wire. To find the current in some wire we need just to 
calculate the difference between the number of particles incoming through this wire and the 
number of particles exiting the central system through the same wire. We do not need to 
know what happened with an electron inside the central system. We only need to know 
the quantum- mechanical scattering amplitudes for an electron to be transmitted/reflected 
through/from the central system. The advantage of the scattering approach is the simplicity 
and the physical transparency of expressions written in terms of scattering amplitudes. We 
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stress that the scattering approach does not aim to calculate the single particle scattering 
amplitudes. This approach tells us how to calculate the transport properties of a mesoscopic 
structure coupled to wires if the scattering amplitudes are known. To calculate the scattering 
amplitudes one can use the Green's functions method. We will give an explicit expression 
for the scattering amplitudes in terms of corresponding Green's functions. Actually the 
combining Green's functions - scattering approach is one of the most powerful and practical 
approaches for transport phenomena in mesoscopic structures. 

1. General formalism 

The scattering approach to transport phenomena in small phase-coherent samples con- 
nected to macroscopic reservoirs was introduced and developed by Landauer and Biittiker 
(Landauer 1957, 1970, 1975, Biittiker 1990, 1992, 1993). 

Within this formalism we consider electrons only in the one-dimensional wires connecting 
the central system to macroscopic reservoirs. It is convenient to introduce separate operators 
aa{s) for incoming and ha{s) for scattered electrons with energy e. 

Then the current Ja{t) flowing into wire a to the central system is the following (Biittiker 
1992): 

Ja{t) = -^J dede'e''^' {{ai{e)a^{e')) - {hi{e)h^{e'))] . (54) 

Here (...) denotes averaging over equilibrium states of reservoirs. 
Correspondingly the dc current reads: 

Ja = j^J Je{Ue)-fi-^He)}, (55) 

where fa^^\e) — {bl,{e)ba{s)) is the distribution function for electrons exiting the central 
system through the wire a, and /«(£) = {al^{e)aa{s)) is the distribution function for electrons 
incoming through the wire a. This expression tells us that the dc current is the difference per 
unit time between the number of electrons entering and exiting the system. As the reservoir is 
at equilibrium, the distribution function for the incoming electrons is the Fermi distribution 
function. In contrast, the scattered electrons, in general, are non equilibrium particles. To 
calculate the distribution function for the scattered electrons we express the 6-operators in 
terms of a-operators. Since an electron coming from any wire can be scattered into a given 
wire a, then the operators 6^ depend on all the operators for the incoming particles. In the 
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model we consider in this chapter, the number of reservoirs is two, /3 = L, R. Therefore, 



normahzed in such a way that their square define corresponding currents (Biittiker 1992). 
The quantities Sapie) can be viewed as the elements of some matrix, which is called the 



2. Floquet scattering matrix 

If the scatterer is driven by external forces which are periodic in time with period r = 
27r/r2o, then interacting with such a scatterer an electron can gain or loss some energy 
quanta nhflo, n = 0, ±1, .... Therefore, in this case the scattering amplitudes in addition to 
the two wire indexes become dependent on the two energies, one for the incoming and the 
other for the outgoing electrons. Such a scattering matrix is called the Floquet scattering 
matrix Sp (see, e.g., Platero and Aguado 2004). Their elements, Sp^aisi^nj^), are related to 
photon-assisted amplitudes for an electron with energy e entering the scatterer through the 
lead f3 and leaving the scatterer with energy En = £ + nMl^ through the lead a. Now the 
relation between the operators h for outgoing particles and a for incoming particles reads 
(Moskalets and Biittiker 2002a): 



where the sum over n runs over those n for which e„ > £0/3, hence it corresponds to propa- 
gating (i.e. current-carrying) states. We denote the Floquet scattering matrix the subma- 
trix corresponding to transitions between the propagating states only. In the case where 
tiTLo <^ e, the sum in Eq. (1561) runs over all the integers: —00 < n < 00. In what follows 
we assume this to be the case. Note that if the scatterer is stationary, the only term that 
remains non-vanishing is that with n = 0, and the Floquet scattering matrix is reduced to 
the stationary scattering matrix with elements Sais^s) = Sp^aisi^,^)- 

The conservation of the particle current at each scattering event implies that the Floquet 
scattering matrix is a unitary matrix (Moskalets and Biittiker 2002a, 2004): 




scattering amplitudes. These amplitudes are 





(56) 



(3=L,R n 



00 




(57) 
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'S'^,a/3(^, £^n) 'S'i;'^7/3(^m; ^n) — ^mO ^0:7 • (58) 

/3=L,_R n.=-oo 

Using Eq. we calculate the distribution function for electrons scattered into wire a: 

00 

fi=L,R ra=— 00 

This function is not the Fermi distribution function unless the scatterer is stationary and 
all the reservoirs have the same chemical potentials and temperatures. This reflects the 
fact that the particles scattered by the dynamical scatterer (quantum pump) are out of 
equilibrium. 



3. Adiabatic scattering 

If the driving forces change slowly, i^o — ^ 0, they behave as if they were almost constant 
for the electrons propagating through the central system. For this reason, the scattering 
properties of a slowly driven (adiabatic) scatterer are close to those of a stationary one. 
Nevertheless there is an essential difference: in spite of the slowly change of the fields, an 
electron can still absorb or emit one or several energy quanta hQo in its travel through the 
central system. Therefore, although the adiabatic scatterer is characterized by the Floquet 
scattering matrix dependent on two energies, SF{£n,£), it is natural to expect that it could 
be related to the stationary scattering matrix S{e) under these conditions 

The stationary scattering matrix S depends on the electron energy e and some properties 
of the scatterer. To account the latter dependence we introduce the set of parameters, 
{Pi}, i = 1, . . . ,Mp and write S{{pi},e). Under the action of external periodic forces the 
parameters periodically change in time, Pi(t) = pi{t + r). Therefore, the matrix S becomes 
dependent on time, S(t,e) = S{{pi(t)},e) and periodic, S(t,e) = S{t + T,e). The obtained 
matrix is called the frozen scattering matrix. This name means that the matrix ^(to,^^) 
describes the scattering properties of a stationary scatterer whose parameters coincide with 
the parameters of a given scatterer frozen at time t = to- The Fourier coefficient for the 
frozen matrix, 

Sn{e)= / -e^"^°*^(t,£), (60) 
can be related to the Floquet scattering matrix. 
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At low driving frequencies, Jlo 0, one can expand the elements of the Floquet scattering 
matrix in powers of Qo. Up to the first order in Qq we have (Moskalets and Biittiker 2004): 

SF{en,e) = S4e) + ""f' ^^^^'^ + hQoAn{e) + 0(^1) . (61) 

Here An is the Fourier transform for a matrix A{t, e), which formally encloses corrections that 
can not be related to the frozen scattering matrix and has to be calculated independently, 
see (Moskalets and Biittiker 2005) for some examples. Note that in the above equation the 
frozen scattering matrix and the matrix A should be kept as energy-independent within a 
scale of order hQ^. 

The unitarity of the Floquet scattering matrix puts some constraint on the matrix A. 
Substituting Eq. ( 16T]) into Eq. ( 157|) and taking into account that the stationary (frozen) 
scattering matrix is unitary we get the following relation: 

.^5t^ + At5}^-(^ — --^-j. (62) 

The advantage of the adiabatic ansatz, Eq. (1611) . is that the matrices S and A depend only 
on one energy and thus have a much smaller number of elements than the Floquet scattering 
matrix. In addition, the adiabatic ansatz allows us to draw some conclusions concerning the 
physical properties of slowly driven systems, in particular, concerning the generated heat 
flows. 



C. Floquet scattering matrix versus Green's function 

There exists a simple relation between the Floquet scattering matrix elements and the 
Fourier coefficients for the Green's function (Arrachea and Moskalets 2006): 

where the Floquet component of the Fourier transformed Green's function Q{n,uj) was in- 
troduced in Eq. 0461) . The equation 0631) is a generalization to periodically driven systems of 
a formula proposed by Fisher and Lee (Fisher and Lee 1981) for stationary systems. This 
relation is based in the fact that the unitary property (1571) and (1581) which is fundamental to 
prove the conservation of the charge within the scattering matrix formalism can be proved 
from identities between the Green's functions, see Eq. (H71) through the relation (!63|l . We do 
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not present in this chapter further details on those proofs. Instead, in the next subsection 
we exphcitly show that both formahsms lead to expressions for the currents through the 
contacts that are equivalent provided that the above relation holds. 



D. Final expressions for the dc currents and powers 

1. Particle currents and the conservation of the charge 

We begin with the expression for the dc particle currents within the Green's function 
formalism. In the subsection fllll A ip we have expressed instantaneous values of observables 



in terms of lesser Green's functions. Now, we use those expressions to evaluate the dc 
components of the currents defined in subsection ( III CI) . In particular, for the charge currents 
ffTOj) and (dH) we have: 

Ji = ^ c?tRe[G<i/t,t)], 

= dtRe[Gl,^{t,t)]. (64) 

Using the representation fj46l) in fl34l) and fl35l) and substituting in the above expressions 
casts for the charge currents within C: 

— 2ew °° f^°° 

Jl = — 5Z X] / dufa{fkj)Ta{u)lm[gi+i^io,{n,uj)gii^{n,uj)], (65) 



a=L,R n=—oo ' 

and through the contacts 



2e\wr 



2 /•+00 



Ja = J dujRe{ifaihuj)gia,laiO,^)Paiuj) 



-co 



Yl Yl 5Z'^/^(^)l^'"''/5(?^,^)|%(^)c/L')fcaK 
f3=L,Rn=—oo ka 



hoo 



duj< fa{fkj)Ta{uj)2lm[gia,ia{0,uj)] 

+ 00 s 

- X] X] fl3{hw)Ta{uJn)\gia,l(s{^^^)\'^^l3{^)\- (66) 
0=L,Rn=-oo 

In the above equations we have used the definitions of the density of states ([3]) and the 
functions (l37Il . Going from the first to the second identity, we have also used the property 
^^[dk'akai^)] ~ ^^^[dk'akai^)] ~ ~Pa{^)/'^, which Can be easily derived just evaluating the 
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Fourier transforms in (136|) . From the identity (H7|) , this current can also be expressed in the 
more compact and symmetric form: 



g ^ r+oo 

Ja = ^ ^ / duj[fa{hWn) - ff3{hhj)]Ta{uJn)\Gla,l(3{'>T'^^)\'^'^pi^) ■ (67) 

l3=L,Rn=-oo-^ 

Within the Floquet scattering matrix approach, we proceed as follows. Substituting 
Eq. (!59|) into Eq. (155|) we get the current in terms of the Floquet scattering matrix elements: 

Ja = T I d<^{ fa{e) - y ] y] //?(£„) \SF,ap{e,en)\^ } . (68) 



h 



/+00 ( go 1 

I P=L,Rn=-oa ) 



An equivalent expression is obtained if we make a shift — > e (under the integration over 
energy) and an inversion n —n (under the corresponding sum) in the term containing 
fp^Sn)- The results is: 



/+00 ( CO "I 

[ 8=L,Rn=-oo J 



J- = lJ de |/,(e) - > : > : Ue) |^i.,„^(e„, e)|^ } . (69) 

Finally, we can write this equation in an alternative way as follows. We multiply the term 
fa{e) in Eq. by the left hand side of the identity, {Sajsie, £n)f = 1, following from 

the unitarity condition Eq. fISS]) . change e„ — > e and n —>■ —n in the resulting expression, 
and find: 

^"=h^ 5Z / de[Ue^)-fp{e)]\SF,^p{er„e)\\ (70) 

l3=L,Rn=-oo-^ 

It is important to note that Eq. ( |68l) coincides with ( |66l) . while ( 1701) coincides with ( 1671) if 
we apply the relation (163|) between the Floquet scattering matrix and the Green's function. 

Another feature worth of being mentioned is the fact that from the expressions (|67|) and 
(ffOl) it can be proved the conservation of the charge, which implies: 

Y.X = 0. (71) 

a=L,R 

We recall that the Jq, was defined as the current exiting the reservoir, for this reason current 
conservation implies that it has different signs at the two reservoirs. A final issue that 
becomes apparent from Eqs. ( 1671) and ( 170i) . is the fact that for slow driving, Qq 0, only 
electrons near the Fermi energy e ^ fi will be excited and hence will contribute to the 
generated current, in agreement with our intuition. 
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2. Particle currents within the adiabatic approximation 



In the subsection (IIIIB 3P we have introduced an approximation for the low driving hmit of 



the full Floquet scattering matrix that depends on the frozen scattering matrix and a matrix 
A. In this section we present the expression for the current in terms of that approximation. 

We have mentioned that the unitary condition imposes a constraint to the matrix A. 
Another more specific constraint follows from the conservation of a charge current expressed 
directly in terms of S and A matrices. To derive it we calculate the dc pumped current 
Ja up to Ql terms for all reservoirs at the same temperature and chemical potential, i.e. 
fa = /ojVo;. Since in the adiabatic case under consideration Qq 0, then at any finite 
temperature it is ksT ^ hflo, and we can expand foi^) — foi^n) ~ —{dfo/de)nhflo — 
{d'^ fo/de'^){nhQo)^ /2. Substituting this expansion and Eq. fl6Tl) into Eq. flTOj) and performing 
the inverse Fourier transformation we calculate the charge current as a sum of linear (upper 
index "(1)" ) and quadratic (upper index "(2)" ) in driving frequency contributions, Ja = 

+ Ja^ + oinl),with 

jf..±r^f.mr^i,j2n,At,.f^] . (73) 



27r i_oo \ de J Jq T \ dt 
The linear behavior of the current as a function of the frequency was calculated by 
Brouwer (Brouwer 1998) using the scattering approach to low-frequency ac transport in 
mesoscopic systems developed by Biittiker et al. (Biittiker et al. 1994). The conservation 
of this current, Yla "^q ^ ~ demonstrated by Avron et al. (Avron et al. 2004) on the 

base of the Birman-Krein relation, (iln(detS') = —TT{SdS'^) (where det(X) and Tr(X) are 
the determinant and the trace of a matrix X, respectively), applied to the frozen matrix 
which is unitary. 

The conservation of the current up to the second order in frequency, =0, leads 

to the constraint for the matrix A we are looking for: 

Equations (16T!) and (162|) show us that the expansion in powers of Qq actually is an 
expansion in powers of /ifio/fe, where 6e is an energy scale characteristic for the stationary 
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scattering matrix. The energy 5e relates to the inverse time spent by an electron with energy 
e inside the scattering region (the dwell time). Therefore, one can say that the adiabatic 
expansion, Eq. (1611) is valid if the period of external forces is large compared with the dwell 
time. It is important to stress that this definition of "adiabaticity" is different from that 
usually used in quantum mechanics one which requires the excitation quantum Ml^ to be 
small compared with the level spacing. 

3. Particle currents within perturbation theory 

In order to gain physical intuition on the behavior of the dc charge current, let us consider 
the weak driving regime (low Vq) and let us evaluate (!67|) with the perturbative solution 
of the Green's function we have presented in fl52l) . We assume that both reservoirs are at 
temperature Tq, = 0. Substituting ( l52l) into ( |671) . we get: 

— e /""^ 

Jo. = -f^^ ^ I du[fa{f^k)~fi3{h^Wa{uJk)\Gi^,ii3{k,uj)\^Tp{uj). (75) 

In the same spirit as in the adiabatic approximation, let us consider that the driving is slow, 
i.e. Qq 0, and let us expand the integrand of the above equation up to the first order in 
Qq. Replacing the Floquet components evaluated up to second order in perturbation theory 
( l53l) we get: 

9 T/2r) ^ 

Ja = ^-p E E r«(/i)r,(/x)sin(5,-5,0Gl!;^,(/^)G^l'i?,(/^) 

j,j'=l (3=L,R 

(76) 

where j,j' runs over the M pumping potentials. Thus, even without specifying the geomet- 
rical details on the structure, which are contained in Eq. (1761) provides us a valuable 
piece of information. As a first point it tells us that at low driving the leading contribution 
to the dc particle current is oc VqQo A second important point is that with local time- 
dependent potentials, as we are considering in our model, we need at least two of these 
potentials operating with a phase lag in order to have a non-vanishing value for this lowest 
order contribution. 
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4- Energy and heat currents 



We can follow a similar procedure as in the previous subsection to derive the dc energy 
and heat currents. In terms of Green's functions we start writing the dc energy currents 
(IT^ and ^ as follows: 

—E 2w r 

Ji = n^J^ dt {Re[G<,,^i{t,t)]w -Re[G<,,^i^,{t,t)]ei+i{t)} , 

= ^J2rdt e,^ Re[G< t)] , (77) 

The energy current within C is: 

a=L,Rn=-oo 

-ei+iIm[^;+2,/„(n,a;)^;+i ,„(n,a;)]} , (78) 

where we have assumed that the position / + 1 does not coincide with a pumping center, 
while for the energy current through the contact we get: 

I |2 p-\-oo ^ 

- ^ ^ fl3{hhj)5{ujn- eka/K)\Qla,l(3{n,Uj)\^Tp{uj)\ 
l3=L,Rn=-oo ^ 

= — duj{ujfa{hi^)Ta{uj)2lm[gia,ia{0,uj)] 

- ^ ^ u;n//3(/ia;)ro(cj)|^/a,i^(ra,u;)pr^(cj) L (79) 

f3=L,Rn=-oo ^ 

which can also be written in the symmetric form: 

= ^ ^ ^ / [/a(/ia;„) - //3(/ia;)]rQ(u;„)|^/a,;^(ra,u;)|^r/3(u;) . (80) 

We now go back to our heuristic argument introduced in Section IIIDI to define the heat 
current. The above equation shows that for low driving, even for reservoirs at T = and 
very weak driving, such that Qq 0, there is a finite energy flow J„ oc fiJa, with given 



in (1671) . This energy is transported by the currents from one reservoir to the other one, thus 
having a convective character and should be subtracted to get a heat flow. 
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To calculate the heat flow we multiply (166!) by /i/e and subtract it to (1791) : 

^ {uJn - UJF)fl3{^)Ta{uJn)\gia,ll3{n,Uj)\'^Ti3{uj)[ , (81) 
l3=L,Rn=-oo ^ 

where hjjp = /x. Equivalently, from (IHOl) and ( |671) . we can write the heat current flowing 
through the contact as follows: 

^ +00 „+oo 

= — ^ ^ / - c^ir)[/a(;ia;„) - /^(;),a;)]rQ(u;„)|^;Q,z/3(n,u;)pr/3(u;) . 

(82) 



f3=L,R n=—oo ' 



Within the scattering matrix approach, one can also calculate the heat current by 
analogy to the charge current, Eq. (1681) . As we already mentioned, Eq. (!68|) contains the 
difference of number of electrons with energy e entering and leaving the scatterer through 
the same wire. Each of these electrons has an energy e. Therefore, to calculate the heat 
current we multiply the integrand in Eq. ( l68i) by {e — fi), drop an electron charge e, and get: 



I r+oo r oo ^ 

I l3=L,Rn=-oo ) 



(83) 



This equation is equivalent to (ISTl) through the relation (163|) . Next we make shifts En ^ £ 
and n — ^ —n in the term containing fi3{en) and flnally obtain: 

1 r+oo r oo ^ 

Ja=-fl dsUe- ^)Ms) - J2 E - \SF,api^n,e)f \ , (84) 

I l3=L,Rn=-oo ) 

We multiply the term containing /^(e) by the identity 1 = 'Yl,pTln \SF,af3{£^£n)\'^ ■ Then we 
make shifts e and n —n in this term and flnally get: 

CO P + OO 

Ja =T Yl Yl / de{en-fi) [fa{£n) " Ms)] \S F^^n, £)\^ , (85) 
f3=L,Rn=-oo^-'^ 

which, because of ( l63l) . is equivalent to (l82l) . 



5. Mean power developed by the fields 

The dc power ( fT7|l done by the ac flelds reads: 

P, = ^[d/jflG<(t,t). (86) 
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In terms of the representation (H6|) it results: 

fen pT/O r+oo 

This expression does not have a counterpart in terms of the Floquet scattering matrix. 
This is because the evaluation of this quantity depends on the microscopic details included 
explicitly in the Hamiltonian. In fact, notice that the formula (IMj) relates the scattering 
matrix only with the Green's function with the coordinates of the central system, la, 1(3, 
that intervene in the contacts. For the same reason, we have shown in the previous section 
equivalent expression within both formalisms only for the currents through the contacts and 
not for the currents within C. Nevertheless the total power developed by all the fields can 
be also calculated within the scattering matrix formalism, see Eqs. ( l88|) and (!97l) . 



E. Technical summary 

To close this section we present in Figure H] a diagram with the summary of the procedure 
to evaluate the different physical quantities we need to discuss the transport behavior of a 
quantum pump, the alternatives and the possible approximations. 



IV. RESULTS AND CRITICAL DISCUSSION 



In this section we apply the concepts and techniques introduced in the previous section 
to analyze the conservation of the energy and the different mechanisms of heat transport 
that we can identify in our quantum pump. On the basis of our previous definitions we can 
show the existence of three generic effects due to a dynamical scatterer. At any segment 
of the system, it is possible to verify the conservation laws introduced in section IIIBI by 
numerically solving the Dyson equation for the retarded Green's functions, evaluating the 
relevant expectation values of observables following the indications of the diagram of Figure 
m In what follows we present analytical results based on the perturbative solution of the 
Green's function and the adiabatic approximation for the scattering matrix. Without the 
explicit evaluation of the functions G^'^^u), which depend only on the geometric statical 
properties of the system, this procedure allow us to analyze the physical properties of our 
system within the weak driving regime. 
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Green's functions 



Scattering matrix 



Eq.{44)for 
Eqs.(37)for 2^(6^) 



Eq.{45]for G^'it^O)) 

Eq.{46)for ^{k,(0) 

Perturbative solution 
Uqs.{52-53) . 



Eq.(63) 



Eq.{65) for 


J, 




Eqs.(66-67) for 






Eq. (77) for 






Eqs.(81-82) for 


a 




^Eq.(87] for 




) 



Eq.(53). 



Eq.(56) for the Floquet 
scattering matrix. 

Adiabatsc approximatiorr 
Eq.{61). 





r 




Eqs.(68-70] for 




Eqs.{83-85] for 


a 


1. 


J 



FIG. 4: Diagram summarizing the possible steps to be followed in order to evaluate particle, energy 
and heat currents as well as the power developed by the fields by recourse to the two formalisms 
presented in this chapter. 

To make the effects clearer we consider the case when the two reservoirs have, not only 
the same electrochemical potential /Xq = /x, but also the same temperature = T ^ 
faihu) = fo{huj),\/a. 



A. Heating of the reservoirs by the quantum pump 

The first effect that takes place in our quantum engine is the heating of the reservoirs (see, 
e.g., Avron et al. 2001, Moskalets and Biittiker 2002, Wang and Wang 2002, Avron et al. 
2004). Unlike the charge current, J^, the sum of heat currents in all the wires, J^^j = J^, 
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Vg cos(Qpf + ^j) 




T=0 



FIG. 5: Scheme of the working regime of the quantum pump when the two reservoirs are at 
temperature T = 0. All the power developed by the external fields is dissipated in the form of heat 
that is absorbed by the left and right reservoirs. 

is non zero. According to the conservation of the energy expressed in Eq. ( 1191) . the definition 
of the heat current (120]) and the conservation of the charge fl71l) . it is clear that the total 
power developed by the fields is equal to the total heat current that enters the reservoirs: 



For reservoirs at temperature T = our intuition suggests us that the total power developed 
by the fields is fully transformed into heat which flows into the reservoirs (see Figure [5]). In 
what follows we analyze the behavior of this flow as a function of the pumping parameters 
within the low driving regime. To this end, we follow an analogous procedure as in subsection 
nil D 3\ and we use perturbation theory to evaluate the powers and heat flows at the contacts. 

1. Heat current at weak driving, T = 

Following exactly the same lines as those presented in the derivation of (1751) . we start 
from ([H2D, we substitute the perturbative solution of the Green's function (15^ and (155]) and 
expand in Taylor series the resulting expression up to the lowest non- vanishing order in Qq. 



M 




(88) 



1=1 



a 
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The result is: 



j,j'=l P=L,R 



The total heat flowing through the contacts reads: 

M 

^^r ■ 

TT 



^ ^ E cos(5, - 5,)V^MT,[u,) 



a=L,R 



j,j'=l 0=L,R 
y COs{6j - 6f)\pij,lf{ujF)\ , 

IT ' ■* 



(89) 



(90) 



where we have used the identity between equilibrium Green's functions and the deflnition 
of the matrix presented in fl5ip . Thus, at T = and weak driving, there is a net heat flow 



oc V^^^^o i^^o reservoirs. 



2. Mean power at weak driving, T = 

We now follow a similar procedure to evaluate the mean power developed by the j-th 
force. Substituting the perturbative solution fl52]) in fISTj) . and keeping terms that contribute 
at 0{V^) we get: 

p. ^ E E / dujfa{huj)Ta{uj) 

xlm {e-^^^[^,,-,„(0,a;)^*.,Jl,a;) + ,„(-!, o;)^*-, JO, a;)]} . (91) 

Then, replacing (l53il we derive an equation with several terms which can be collected as 
follows: 

M 



j'=i 



with 



= ^^o(^^o) ^ d^/o(;«a;)Im{7jj/(w)7j.^.,(w)} 



oo 
oo 



Ag, = ^^M^ f °°d^/o(MRe{7,,,'(^)7iv(^)} , (93) 

^ J —oo 
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being 



a=L,R 



la 



yO,R 



(94) 



3. Conservation of the energy 

The second term of fl92l) vanishes when we perform a summation over all the fields, since 

(2) 

\ - -, is symmetric under a permutation j ^ j' while sin(5j — 6j) is antisymmetric under this 
operation. Thus, the only term contributing to the sum over all the powers is the first one, 
which for low Qq results: 




xRe I [fo{hw - hno)pij,ij'{^^ - ^o) - /o(^ + hno)pij,ij'{^^ + ^o)][G'°jJ,(a;)]*| 

(95) 

Performing the sum over j in (l92l) and using \pi^if{uj)\'^ = [^^^^'(ci;)^ + \Gi'j{uj)\'^ — 
2'Re[Gi^ir{uj)Gir^i{uj)], we can verify the fundamental law of the conservation of the energy 



evo) 



TT 



-Plj,lj'{^F) G^IuA^f) 



B. Energy exchange between external forces 

The evaluation of the coefficient A^- j, at weak driving can be carried out following exactly 
the same steps as with A^^j,. The result is 

\2 /•+00 



(2) 2/lf2o(eVo) f °° J J lr^o,R i \r^o,R , \1 /n«^ 
ij' / ^^^"^ ' (96) 

" J —oo 



^3,3' 



i.e. this contribution is oc VLq, and therefore dominates the behavior of Pj at weak driving. 
Interestingly, this contribution does not exist in a configuration with a single ac field, while 
it can have different signs at different fields in a configuration with several pumping centers. 

Therefore, we present the second general effect taking place in quantum engines: One 
external force can perform work directly against another external force with a negligible 
amount of energy being dissipated into the reservoirs. This remarkable mechanism opens 
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Vq cos(D.Qt -\- Si) 






FIG. 6: Scheme of the working regime of the quantum pump when the two reservoirs are at 
temperature T = and low driving: Vq and f^o small. The dissipated energy flowing into the 
reservoirs is low, while it is possible that part of the work done by one of the ac fields is coherently 
transferred to the the other one, which receives the ensuing energy. 



the possibility of the coherent energy transfer between pumping centers as indicated in 
Figure [61 



C. Directed heat transport at finite temperature 

To show that the dynamical scatterer can induce a directed heat transfer between the 
reservoirs we, first, calculate the total generated heat J^t = J2a "^S- Summing up Eq. (18^ 
over a we find (for = /o,Va;): 

Q °° POO 

J?ot = -i^Y. E E / deUe)n\S,^^p{e^,e)\\ (97) 

a=L,R I3=L,R n=-oo 

The part of the total generated heat which flows into wire a, jf^t = J2a '^a,geni can be 
defined as follows: 

4^,en = -^E E / defoie)n\SF,^pien.e)\\ (98) 

D ^ "'-00 



f3=L,R n=— 00 ' 

The remaining part of the heat flowing into wire a, J^^pump — ~ J a, gem i^: 



^ r+00 r 00 ^ 

jS,pump =1 de{e - /i) /o(£) <^ Y \SF,ap{en, £) |' " U 

Kf3=L,Rn=-oQ ) 



(99) 
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FIG. 7: Scheme of the working regime of the quantum pump when the two reservoirs are at a finite 
temperature T. There is a net pumping of heat from the one reservoir to the other. The quantum 
pump, thus works as a refrigerator. 

Using the unitarity condition for the Floquet scattering matrix, Eq. fl57|) . one can easily show 
that the part of the heat current J2,pump satisfies the conservation law similar to the one for 
the charge dc current, Eq. (17T]) : 

E Jivump = ^- (100) 

a=L,R 

This means that J ^ pump is transported from one reservoir to another one with the help of 
a dynamical scatterer. By analogy with the corresponding charge current we identify this 
portion of the total heat as a pumped heat (hence the lower index "pwrnp"). This is the 
third general effect we identified in our quantum engine: The dynamical scatterer induces a 
directed heat transport between the reservoirs (see, e.g., Humphrey et al. 2001. Segal and 
Nitzan 2006, Arrachea et al. 2007, Rey et al. 2007, Martinez and Hu 2007). 

If the pumped heat is, for instance, negative in the L wire, J^^pump < 0, then it is 
necessarily positive in another wire, J^^pump > 0. If the absolute value of this heat is larger 
than the one of the generated component ^g^, then the whole heat flowing into the R wire 
is positive, i.e. directed from the reservoir to the central system, = + J'^p^mp > 0. 

In this case the reservoir R will be cooled while L will be heated. 

The splitting of into Ja.gen Ja,pump helped us to show that can be positive. 

38 



Strictly speaking, such a splitting is not unique and only the whole heat current has a 
direct physical meaning. However at slow driving, one can support such a decomposition of 
into the generated and the pumped heat by additional physical arguments as follows. 



1. Adiabatic heat currents 

The expansion ( I6T1) allows us to calculate the heat flow with an accuracy of 0{Q'^). To 
show it explicitly we rewrite slightly Eq. (18^ (with /„ = /o, Va). We assume fc^T ^ hflo 
and expand the difference of Fermi distribution functions in fl5^ in powers of Qq, use Eq. flM]) 
and find from Eq. §E) the heat current, = J^'^^^ + jS'^^^ + O (fig), where 



7Q,(2)__A (-^\ r dt [ dS{t,e)dS^{t,e 

" " 4vr y_ dejJo r[ dt dt 



~ I ^de{e-fi) r-lm[2QoA{t,e)^^^^^] . (102) 




2vr J_^ V de 

Next we split the heat current into the generated heat and the pumped heat as follows, 

jQ ^ jQ _|_ jQ ™4-v, 

a a,gen ' " a, pump i "^'■^'■^ 



aa 



4«„ = -^/;j<ie(-.)(-f)/'flm([s + 2RM]^) . (104) 

Notice that these equations also remain valid at ultralow temperatures, ksT <^ hQo, which 
can be verified by direct calculations taking into account the energy-independence of the 
matrices S and A over a scale of order Qq, i.e. over the region of the thermal widening of 
the edge of the Fermi distribution function. 

The above given splitting is justified by the following observations, (i) The quantity J^^en 
is negative in each wire a as it should be for the heat generated by the scatterer and flowing 
into the reservoirs, (ii) At zero temperature the pumped heat vanishes identically, Ja,pump — 
0, since it is impossible to take heat out of the system kept at zero temperature. To prove 
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the first observation we show that the integrand in Eq. fllOSp is positive. To this end we use 
the Fourier transformation and get, 1/r dt{dS /dt dSy dt)aa = 1^/3 ^^l*^"/?,"!^ > 0- 
The second observation follows from the fact that at zero temperature it is (e—^) dfo/de = 0, 
hence the equation fll04p vanishes. Note that the conservation of the pumped heat current, 
J2a Ja,pump — 0; dircctly follows from the conservation of charge currents, Eqs.([72]) and (1751) . 
which implies 1/r ^ dtlmTi S + 2hVLQA dS^dt = 0. 

From Eq. fllOSp it follows that the adiabatic scatterer heats the reservoirs with a rate 
proportional to hQ^ (Avron et al. 2001). In contrast, the pumped heat, Eq. fll04p . is rather 
proportional to ksTVLQ. At sizable temperatures, ksT ^ hflo, the amount of pumped heat 
can exceed the generated heat, \Ja,pump\/ Ja,gen ~ ksT/^hQQ) ^ 1. Therefore, if in the wire 
a we have Ja,pump > 0, then the reservoir a will be cooled (see Figure [7]). This mechanism 
opens the possibility of using quantum pumps as refrigerators. 



V. SUMMARY 



In this chapter we have introduced the basic concepts to analyze at the microscopic level 
the energy transport in quantum systems driven by harmonically time-dependent fields. We 
have introduced a simple microscopic model for a quantum pump, which consists in a finite 
structure connected to two macroscopic reservoirs, with ac local fields that oscillate in time 
with the same frequency and a phase lag. We have analyzed the fundamental conservation 
laws for the charge and the energy and we have defined the basic concepts to study the trans- 
port behavior in these systems: charge currents, energy currents, heat currents and powers 
developed by the fields. We have reviewed two complementary techniques to calculate the 
currents and the powers: the non-equilibrium Green's function formalism for harmonically 
time-dependent Hamiltonians and the scattering formalism for periodically driven meso- 
scopic systems. We have shown that the two approaches are equivalent for the evaluation 
of the charge and heat currents through the contacts between the driven system and the 
reservoirs. We have also introduced two approximations: the adiabatic approximation to 
the Floquet scattering matrix and a perturbative solution of the Dyson's equations for the 
Green's functions valid within the weak driving regime. Both techniques are important to 
draw conclusions on general features of the transport behavior without the explicit evalua- 
tion of the Green's functions or the scattering matrix elements. Such conclusions are, thus. 
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generic and do not depend on the geometrical details of the driven structure. A summary 
of the technical details, including the main equations and the alternative routes to evaluate 
them exactly or in an approximate way is given in a diagram at the end of section III. 

Finally, in section IV we have applied the concepts and tools we have introduced in 
the previous sections in order to discuss three important mechanisms of energy transport 
in quantum pumps. The first one is the fact that the total work done by all the local 
fields is dissipated in the form of heat that flows to the reservoirs. This effect is rather 
expected. In any case, we have exploited our theoretical techniques at weak driving to 
evaluate term by term powers and heat currents and explicitly verify the conservation of 
the energy. To unveil a fundamental law is always a beautiful result in theoretical Physics 
and an important support for the power of a theoretical tool. In addition wc have shown 
that other two less expected and subtle transport mechanisms can take place: the coherent 
transport of energy allowing for regimes where some of the forces make work, while other 
receive work. This interesting mechanism could be exploited, for instance, to couple two 
quantum pumps in a combined engine. The final remarkable mechanism is the pumping 
of heat at finite temperature and weak driving, allowing for the operation of the quantum 
pump as a refrigerator which extracts heat from a reservoir and injects heat in the other 
one. 

VI. FUTURE PERSPECTIVE 

The different operational regimes that we have identified in the quantum pumps have 
several important outcomes. On the theoretical side there are several lines to further analyze. 
A first issue to explore is the role of the geometrical details of the structure, in order to 
identify the optimal architecture to enhance each mechanism and improve the efficiency of 
the quantum engine. Another important ingredient is the investigation of the role of many- 
body interactions. In particular, the electron-electron and the electron-phonon interactions. 
On the experimental side it would be very interesting the design of an experimental setup 
to implement these effects. In this sense, it is very promising that quantum refrigeration 
has been already experimentally explored in mesoscopic structures with superconducting 
elements under ac driving [Giazzoto 2006]. 
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